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Abstract 

We argue that in the small x processes, in the black disc QCD regime (BDR) a very 
forward parton propagating through the nuclear matter should loose a significant and 
increasing with energy and atomic number fraction of its initial energy as a result of 
dominance of inelastic interactions, causality and energy-momentum conservation. We 
evaluate these energy losses and find them to lead to the significant suppression of the 
forward jet production in the central NA collisions at collider energies with a moderate 
suppression of recoiling jet at central rapidities. We confront our expectations with 
the recent RHIC data of the STAR collaboration on the probability, P, for emission of 
at least one fast hadron at a central rapidity in association with production of a very 
forward high pt neutral pion in pp and dAu collisions. We calculate the A-dependence 
of P, and find that the data imply a strong suppression of leading pion production at 
central impact parameters. We also conclude that production of recoil jets in the hard 
subprocess is not suppressed providing further evidence for the dominance of peripheral 
collisions. Both features of the data are consistent with the onset of BDR. We suggest 
new phenomena and new observables to investigate BDR at RHIC and LHC. 



1 Introduction 

It is well understood now that one of distinctive properties of hard processes in pQCD is 
the fast increase with energy of cross sections of hard inelastic processes and their signifi- 
cant value. Thus the interactions of the partons produced in the sufficiently small x hard 
processes should be highly inelastic. Dominance of inelastic processes leads to the specific 
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pattern of energy losses for a parton propagating through the nuclear medium which is 
the main subject of this paper. Really in the elastic rescatterings which dominate in the 
large x processes energetic parton looses a finite energy [1] while propagating a distance 
L: AE 0.02GeV L'^ / Fm^ . The analysis of Ref. P] of the irA Drell-Yan pair production 
indicate that the data are consistent with the rate of energy loss by a quark of Ref. [Tj and 
correspond to a energy loss < AGeV for quark of energy ~ 200 GeV propagating through 
a center of a heavy nucleus. In contrast in the deep inelastic processes for example DIS 
off a proton the fraction of initial photon energy lost by incident parton is ~ 10% within 
DGLAP approximations, cf. discussion in section 2. Numbers are probably similar within 
the NLO BFKL approximation corresponding to the rapidity interval between the leading 
particle and next rung in the ladder of about two. (It is equal to zero within the LO BFKL 
approximation which systematically neglects the loss of energy by energetic particles.) 

In the black disc regime the contrast between the different patterns of energy losses 
becomes dramatic. A parton with energy E propagating sufficiently large distance L through 
the nuclear media should loose energy: 



with c ~ 0.1 in small x processes. This energy loss exceeds by orders of magnitude the losses 
in the large x regime. 

Another subtle effect characteristic for a quantum field theory has been found long be- 
fore the advent of QCD: eikonal interactions of energetic particle are cancelled out as the 
consequence of causahty [3l H] . This cancellation including additional suppression of eikonal 
diagrams due to energy-momentum conservation is valid for the exchanges by pQCD ladders 
with vacuum quantum numbers in the crossed channel [5]. The cancellation of the con- 
tribution of eikonal diagrams has been demonstrated also for the exchanges by color octet 
ladders as the consequence of bootstrip condition for the reggeized gluon [6]. Thus suffi- 
ciently energetic parton may experience only one inelastic collision. To produce n inelastic 
collisions wave function of energetic parton should develop component containing at least 
n constituents [5J. This effect leads to the additional depletion of the spectrum of leading 
partons in the kinematics close to BDR where inelastic interactions of the energetic parton 
is important part of unitarization of amplitudes of hard processes. 

Since the number of inelastic collisions is controlled by the number of scattering centres 
at given impact parameter the effect of the suppression of the yield of leading partons should 
be largest at the central impact parameters. We evaluate energy losses of leading parton in 
small X regime of QCD and show that blackening of pQCD interaction leads to dominance 
of peripheral collisions in the production of the leading hadrons/jets in high energy hadron 
- nucleus interactions and to a significant, increasing with energy and atomic number loss 
of finite fraction of leading parton energy in the central collisions. Inclusive cross section is 
cx A^/^ deep in the BDR region with suppression of the recoil jets depending on x of jet. 
One of characteristic features of BDR regime is that there is no suppression of recoil jet in 
the peripheral collisions. At moderately small x which are reached at RHIC, suppression 
of recoil jet should depend on its rapidity and be maximal if both jets carry a significant 
fraction of the projectile energy. We will show that this prediction is supported by the recent 
RHIC data on leading hadron production in dA collisions. 



AE = cE{L/3Fm) 
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It is instructive to compare the kinematics of partons involved in the production of leading 
hadrons at RHIC with that for small x phenomena at HERA. Taking for example the STAR 
highest rapidity (y=4) and {pt) = l-3GeV/c bin [7J we find that xn > 0.7 for the incoming 
parton. Hence, minimal Xg resolved by such a parton are ~ 4:p^/{x]\fS]\fj\f) ~ (2 ^ 3) ■ 10~^. 
This is very close to the kinematics reached at HERA. The analyses of the HERA data within 
the dipole model approximation show that the partial amplitude for the quark interaction 
reaches at HERA strength up to 1/2 of the maximal strength, see review in Ref. [8]. In the 
case of heavy nuclei one gets an enhancement factor ~ 0.5A^^^ so the quark interaction with 
heavy nuclei should be close to BDR for p'f < 1.5GeV'^ and Xprojecuie ~ 0.5. In the LHC 
kinematics BDR will cover much larger pf range, see for example Fig. 17 in Ref. [8]. 

First evidence for suppression of the forward spectra in the deuteron-gold collisions in the 
kinematics rather close to the BDR was reported by the BRAHMS [9j , and further studied 
by PHENIX [10], and STAR[7]. High pt spectra of h~ at 2 < y < 3.2 are suppressed by a 
factor [9] 

dycPpT I dy(P-pT 

which is ~ 0.8 for y = 3.2, pt = 2 GeV/c. Since in the kinematics of the experiment a{pp — *• 
h~ + X)/a{pp + X) < 2, the vr" yield produced by the proton projectile relative to 

that for the deuteron projectile (per nucleon) is substantially smaller. As a result in the 
case of the 7r° production which is produced with equal strength by protons and neutrons 
one expects a bigger suppression. For example RdAuiv — '^■'^iVt = 2GeV/c) ^ 0.55 [llj). 
This suppression factor is significantly larger than expected suppression due to the leading 
twist nuclear shadowing. Suppression was observed in the kinematics where the hadron 
production in pp collisions is in a reasonable agreement with the recent pQCD calculations 
based on the NLO DGLAP approximation [12]. Very recently STAR [7] has reported new 
results for the vr*^ ratios for y ~ 4 and pt < 2.0 GeV. They observed a larger suppression 
- RdAu ~ 1/3, which is consistent with a linear extrapolation of R^au io y = 4 taking into 
account the 2/3 factor due to the isospin effects [TT] . 

The STAR experiment also reported the first observation of the correlations between 
the forward 7r° production with the production of the hadrons at the central rapidities 
l^hl ^ 0.75. Such correlations provide a new information about the mechanism of the 
suppression of the inclusive spectrum. 

The discussion is organized as follows. In section 2 we evaluate energy losses of leading 
partons of the proton propagating through the nuclear medium in the kinematics of the 
onset of the BDR and find them to be > 10% for central proton-heavy ion collisions in the 
RHIC kinematics. In section 3 we discuss expectations of the QCD BDR for the spectra of 
the leading particles. In section 4 to disentagle interplay of soft and hard QCD phenomena 
we evaluate correlations between forward and central hadron production using information 
obtained in the BRAHMS experiment ^ on the dependence of the central multiplicity on 
the number of the wounded nucleons. We calculate the dependence of correlation parameters 
studied by STAR, on the number of wounded nucleons and find that the data require this 
number to be ~ 3, which is significantly smaller than the number of wounded nucleons for 
central impact parameters ~ 13 strongly suggesting dominance of the peripheral collisions 
in 7r° production. We also want to stress that the paper considers the yield of partons with 
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transverse momenta < than that typical for the BDR. At the same time pion production 
with transverse momenta significantly larger than that typical for BDR should be domi- 
nated by the scattering at central impact parameters. For example, Color Glass Condensate 
(CGC) inspired models predict for this case enhancement of production at central impact 
parameters by the factor of ~ A^^^ [13]. In section 5 we perform a detailed analysis of the 
correlation observables within peripheral models of pion production constrained to repro- 
duce the inclusive data. We reproduce the observed values of the correlation parameters and 
find that the suppression of the correlation parameter related to production of recoil jets 
observed by STAR [7] is due to soft interactions and does not indicate suppression of the 
pQCD mechanism of the production of the recoil jets. Thus RHIC data are consistent with 
the pattern expected energy losses in central collisions, cf. Ref. [H]. In section 6 we sug- 
gest several new observables which could allow to diminish model dependence of comparison 
between the hard components of the interaction in pp and dAu cases, quantitative study of 
the suppression on the number of wounded nucleons, which also will provide a probe of the 
color transparency effects as well as effects of large gluon fields. 

2 Energy losses of forward parton in the vicinity of 
black disk regime 

Energy losses for the parton propagation through the nucleus medium are dominated in 
moderate x processes by its elastic rescatterings off the constituents of the media due to the 
Coulomb gluon exchange. Therefore they depend weakly on energy and proportional to the 
square of distance propagated by the parton [1]. 

However, the amplitude with color octet quantum numbers decreases with energy due to 
the gluon reggeization in pQCD as [HI [T3] : 



where P(t) is the gluon Regge trajectory with P{t = 0) < 1. Infrared divergences of P{t) are 
regulated by hadron wave functions. At the same time the amplitude due to exchange by a 
ladder with the vacuum quantum numbers in the crossed channel rapidly grows with energy: 



where X{t = 0) ~ 0.2. (For the simplicity we restrict ourselves here by the phenomenological 
fit to the theoretical formulae and to the HERA data on structure functions of a proton.) 
Hence such amplitudes (modeled at moderately small x as the two gluon exchange ladder) 
fastly exceed single gluon exchange term and at larger energies achieve maximum values 
permitted by probability conservation. 

Thus dominance of elastic collisions breaks down at high energies leading to the regime 
where incoherent processes and incoherent energy losses dominate leading to the loss of finite 
fraction of initial energy of a parton, cf. Ref. [16j. This is the major difference from moderate 
X processes considered in [Ij where coherent energy losses seems to dominate. Consequently, 
single inelastic collision of the parton produced in a hard high energy NN collision off another 
nucleon is described by the imaginary part of the two gluon ladder with the vacuum quantum 



Ag oc a^s^W {i + tan(7r/3(t)/2)) 



(3) 



A oc {t + tan((7r/2)A(t))) 



(4) 
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numbers. By definition, the inelastic cross section is calculable in terms of the probability 
of inelastic interaction, Pineiip) of a parton with a target at a given impact parameter h [17]: 



CTmeZ = / d%Pinelib,S,Q'^) (5) 



Since amei is calculable in QCD above equation helps to calculate Pinei{b, s,Q^) . The 
probability of inelastic interaction of a quark is cf . |19l |8] : 

TT^ A 

Pineiib,x,Q'^) = YOis{kt)-j^xGAix,Q^,b), (6) 

where x ~ ik^/sqN, ~ 4fcj , A ~ 2 (for the gluon case Pmdib) is 9/4 times larger). We use 
gluon density of the nucleus in impact parameter space, Q^, 6) (/ d%G a^x^Q"^ ,b) = 

Ga{x, Q"^)) ■ Above equation for the probability of inelastic interaction is valid only for the 
onset of BDR when Pinei{b, s, Q^) < 1 (which is the unitarity limit for Pinei{b, s, Q^)). 

If Pineiib,x,Q'^) as given by Eq|6] approaches one or exceeds one it means that average 
number of inelastic interactions, N{b) becomes larger than one. Denoting as Gcr{x,Q'^,b) 
for which Pinei{b) reaches one we can evaluate N{b,x,Q'^) as 

iV(6, x, Q2) = Ga{x, Q\ b)/GUx, Q\ b). (7) 

As soon as Pimi becomes close to one, we can easily evaluate lower boundary for the 
energy losses arising from the single inelastic interaction of a parton. This boundary follows 
from the general properties of the parton ladder. Really, the loss of finite fraction of incident 
parton energy -e arises from the processes of parton fragmentation into mass M which does 
not increase with energy. For binary collision = T^^z^y- For the contribution of small 
e < 1/4 

e « k\lM^ (8) 

Here fcj is transverse momentum of incident parton after inelastic collision. The spectrum 
over the masses in the single ladder approximation (NLO DGLAP and BFKL approxima- 
tions) is as follows 

rfa oc / dM^/M\s/M^fe{M^ - Akl), (9) 



where we accounted for the high energy behavior of the two gluon ladder amplitude Eq.(j4]). 
We effectively take into account the energy momentum conservation i.e. NLO effects. Con- 
sequently the average energy loss (for the contribution of relatively small energy losses 
(e < 7 ~ 1/4) where approximation of Eq.([8]) is valid): 

For the realistic case 7 = 1/4, A = 0.2 this calculation gives the fractional energy loss of 
6%. This is lower limit since we neglect here a significant contribution of larger e (it will be 
calculated elsewhere). 

In the kinematics of onset of BDR effective number of inelastic interactions becomes 
significantly larger than 1 so in the evaluation of fractional energy loss one should multiply 
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evaluated above energy loss by the factor: N{b). In addition one should account for the 
phenomenon specific for a quantum field theory in small x regime. The sum of Feynman 
diagrams which leads to eikonal contribution at moderately small x is cancelled out at large 
energies as the consequence of the causahty i.e. analytic properties of amplitudes and their 
decrease with the parton virtuality, (cf. Refs. 0,11] and for the generalization to QCD PlS]) 
and/or energy-momentum conservation, cf. ^[5]. At central impact parameters absorption 
at high energies is due to N{b) > 1 inelastic collisions (interaction with several ladders). 
The energy of initial parton is shared before collisions at least between N constituents in the 
wave function of the incident parton to satisfy causality and energy-momentum conservation. 
This quantum field theory effect which is absent in the framework of eikonal approximation 
can be interpreted as an additional energy loss [8]: 

eA{b) ^ iV(6)ejv. (11) 

Here is the energy lost due to exchange by one ladder - Eq. (ITU]) . Above we do not 
subtract scattering off nucleon since our interest in the paper is in energy losses specific for 
nuclear processes in the regime when interaction with a single nucleon is still far from the 
BDR. If collision energies are far from BDR, the energy losses estimated above should be 
multiplied by small probability of secondary interactions. Inclusion of enhanced "pomeron" 
diagrams will not change our conclusions based on the necessity to account for the energy- 
momentum conservation law. 

Yields of leading hadrons carrying fraction of projectile momentum > Xp are rapidly 
decreasing with x^- as oc (1 — xp)"-. For pion production n ~ 5 6. Obviously for large xp 
average values of x for progenitor parton are even larger, leading to strong amplification of 
the suppression due to the energy losses. The spectrum of leading pions is given in pQCD by 
the convolution of the quark structure function, oc (1 — x)",n ~ 3.5 and the fragmentation 
function oc (1 — 2;)™,m~ 1.5 -j- 2 leading to a very steep dependence on xp, oc (1 — x^?) ^. 
As a result for the STAR kinematics x ~ 0.7 and z ~ 0.8 correspondingly energy losses 
of 10% lead to a suppression roughly by a factor [(0.9 — Xp)/{1 — xp)]^. For xp = 1/2 
this corresponds to suppression by a factor of four. In particular, introducing the energy 
loss of ~ 6% in the NLO calculation of the pion production is sufficient [11] to reproduce 
the suppression observed by BRAHMS [9j. Similar estimate shows that average losses of 
~ 8 -i- 10% reproduce the suppression of the inclusive yield observed by STAR[7]. This 
value is of the same magnitude as the above estimate. Also, Eq. (fTT!) leads to much stronger 
suppression for production at central impact parameters than in peripheral collisions. 

In the kinematics of LHC the same kt{BDR) would be reached at x^v which are smaller 
by a factor srhic/ slhc ~ 10~^, while for the same xn one expects much larger values of 
kt{BDR) (see e.g. Fig. 17 in [8j). Thus in the kinematics of LHC the regime of large energy 
losses should extend to smaller xn- 

There are two effects associated with the interaction of partons in the BDR - one is an 
increase of the transverse momenta of the partons and another is the loss of the fraction of 
the longitudinal momentum [16]. The net result is that distribution of the leading hadrons 
should drop much stronger with xp than in the CGC models [20j where only kt broadening, 
change of the resolution scale and suppression of coalescence of partons in the final state but 
not the absorption and related energy losses were taken into account. At the same time, the 
kt distribution for fixed xp should be broader. Note here that the leading particle yield due 
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to the scattering with kt ^ kBOR is not suppressed and may give a significant contribution 
at smaller kt via fragmentation processes. 

This discussion shows that selection in the final state of the leading hadron {xp > 0.3-=-0.5 
at RHIC) with moderately large kt should strongly enhance the relative contribution of the 
peripheral collisions where BDR effects are much smaller. We will demonstrate below that 
these expectations are consistent with the STAR data. 

At extremely high energies where kinematics of the BDR will be achieved for a broad 
range of the projectile's parton light-cone fractions and virtualities, QCD predicts dominance 
of scattering off the nuclear edge leading to: 

, / , CXAV3, (12) 

for a large enough xn and and a wide range of p*. With increase of incident energy the range 
of for fixed x^ would increase. Also the suppression for a given pt would be extended to 
smaller x^. 



3 Interaction of leading partons with opaque nuclear 
medium 

At high energies leading partons with light cone momentum XN,Pt are formed before nucleus 
and can be considered as plane wave if 

{xNs/mN)il/M^) > 2Ra. (13) 

Here M is the mass of parton pair (and bremstrahlung gluon) produced in the hard collision. 
If sufficiently small x are resolved, the BDR regime would be reached: 

ApI/xns < x{BDR). (14) 

In the BDR interaction at impact parameters h < Ra is strongly absorptive as the 
medium is opaque. As a result, interaction of leading parton lead to a hole of radius Ra 
in the wave function describing incident parton. Correspondingly, propagation of parton at 
large impact parameters leads to elastic scattering - an analogue of the Fraunhofer diffraction 
of light off the black screen. However since the parton belongs to a nucleon, the diffraction 
for impact parameters larger than RA + Vstr (where Vstr is the radius of the strong interaction) 
will lead to the proton in the final state - elastic p A scattering. Only for impact parameters 
Ra + '"'str > b > Ra the parton may survive to emerge in the final state and fragment into 
the leading hadron. Cross section of such diffraction is 27r RAVstr- Another contribution is 
due to the propagation of the parton through the media. This contribution is suppressed 
due to fractional energy losses which increase with the increase of energy, leading to gradual 
decrease of the relative contribution of the inelastic mechanism (see discussion in section 5). 

Thus we predict that in the kinematics when BDR is achieved in pA but not in pN scat- 
tering, the hadron inclusive cross section should be given by the sum of two terms - scattering 
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from the nucleus edge which has the same momentum dependence as the elementary cross 
section and scattering off the opaque media which occurs with large energy losses: 



da{d + A ^ h + X) / dxh<fpt 



c,A^'^ + C2{A)A^'^ 



(15) 



dij{d + p ^ h + X) I dxhd^Pi 



The coefficient c\ is essentially given by the geometry of the nucleus edge - cross section for 
a projectile nucleon to be involved in an inelastic interaction with a single nucleon of the 
target. Coefficient C2(A) includes a factor due to large energy losses and hence it decreases 
with increase of the incident energy for fixed Xh,Pt- Deep in the BDR the factor C2{A) would 
be small enough, so that the periphery term would dominate. 

It is worth to compare outlined pattern of interaction in the BDR with the expectations 
of the CGC models for small x hard processes in the kinematics where transverse momenta 
of partons significantly larger than that characteristic for BDR. These models employ the 
LO BFKL approximation with saturation model [21] used as initial condition of evolution in 
ln{xo/x), see [22] and references therein. In these models the dependence on atomic number 
is hidden in the "saturation scale" and in the blackness of interaction at this scale. In this 
model partons interact with maximal strength at small impact parameters without significant 
loss of energy. Note that leading parton looses significant fraction of incident energy in the 
NLO BFKL approximation but not in LO BFKL [23j. As a result the cross section is 
dominated by the scattering at small impact parameters and depends on A at energies of 
RHIC approximately as A^/^[T3]. Also, the process which dominates in this model at central 
impact parameters is the scattering off the mean field leading (in difference from BDR where 
DGLAP approximation dominates in the peripheral processes in the kinematics of RHIC) to 
events without balancing jets. With increase of jet transverse momenta interaction becomes 
less opaque, leading to a graduate decrease of the probability of inelastic collisions and hence 
to the dominance of the volume term. 

A natural way to distinguish between these possibilities is to study correlations between 
production of forward high pt hadrons and production of hadrons at central rapidities. First 
such study was undertaken by the STAR experiment |7]. 

4 Hadron production in soft nucleon-nucleus interac- 
tions at central rapidities. 

The STAR experiment reported correlations between the leading pion trigger and central 
leading charged hadron production. The procedure picks a midrapidity track with \7]h\ < 0.75 
with the highest pt > 0.5 GeV/c and computes the azimuthal angle difference A0 = — 
for each event. This provides a coincidence probability /(A0). It is fitted as a sum of 
two terms - a background term, B/27T, which is independent of A0 and the correlation term 
5'(A0) which is peaked at A0 = vr. By construction. 



We will argue below that the A-dependence of B and 5* is sensitive to dependence of the 
leading pion production on the centrality of the collision. 




/ f{A(f))dA(f) = B+ S{A(j))dA(j) = B + S <1 



(16) 
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The low pt (soft) particle contribution which is uncorrelated in with the trigger orig- 
inates both from the collisions of the second nucleon of the deuteron with the nucleus and 
from interactions with nucleon involved in the hard collision with several nucleons. This 
contribution should grow with A since the low pt hadron multiplicity for ?/ ~ increases 
with A. To make quantitative estimate of this contribution we will make an approximation 
that the rate of these soft processes is weakly correlated with production of the forward pion 
provided we compare the processes at the same impact parameter. This natural assumption 
is valid in a wide range of models including CGC models. It is consistent also with the 
information provided by STAR on the weak dependence of the central multiplicity on of 
the trigger pion, and lack of long range rapidity correlations for low pt processes which was 
observed in many studies of hadron-hadron collisions. Note that at the LHC energies one 
would have to correct this approximation for the correlation of soft and hard interactions 
in the elementary interactions due to more localized transverse distribution of the valence 
partons, see discussion in Ref. [8]. 

Based on generic geometric considerations one expects that the multiplicity should be 
a function of the number of nucleons on the projectile nucleon impact parameter. Within 
this approximation to estimate effects of soft production on the correlation observables we 
can use information on the impact parameter dependence of the hadron multiplicity which 
is available from several dA RHIC experiments. 

Using the BRAHMS data [9j we find that for the STAR cuts can be roughly described 
by a simple parametrization 

R'=[^y, (17) 

with r ~ 0.2. Here the factor of two in the denominator takes into account that each of the 
nucleons of the deuteron experiences, on average, equal numbers of collisions. For example, 
for an average number of collisions Ncou ~ 7.2, Eq. fll7l) gives = 0.77 while the BRAHMS 
data reports R'' = 0.7 ^ 0.75. 

Note in passing, that the Gribov-Glauber approximation for the hadron - nucleus scat- 
tering combined with AGK cutting rules [23] which neglects energy conservation leads to 
R^ = 1. If one takes into account energy conservation - the split of the energy between Ncou, 
and the increase of the central multiplicity with energy oc one roughly reproduces Eq. 

m- 

First we want to find out what information about centrality of the interactions leading to 
production of the leading pion is contained in the A-dependence of B, the probability that a 
fast hadron within the experimental cuts does not belong to the recoil jet. Obviously, with 
an increase in the number of nucleons in the nucleus involved in the interactions practically 
all events would contain at least one particle in the cuts of STAR leading to B very close 
to one even if the elementary hard interaction is not affected by the nuclear environment. 
Using Eq. (fT7|) we can express B for collisions with n nucleons, Bn through characteristics 
measured for pp collisions. 

Let us denote the probabilities to produce a hadron within the central cuts of STAR due 
to soft and hard interactions in pp collisions hy pb and ps respectively. Since the px cut of 
STAR is rather high (comparable to the momentum of the leading hadron in the recoiling 
jet for the trigger jet with (pt) ~ 1.3 GeV/c), we will assume that in the pp events where 
both soft and hard mechanisms resulted in the production of a hadron (hadrons) within the 
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STAR cuts there is an equal probability for the fastest hadron to belong to either the soft 
or hard component (this is essentially an assumption of a reasonably quick convergence of 
the integrals over pt for prmin = 0.5 GeV/c). [J Within this assumption, the probability 
to produce no fast hadrons is (1 — pb){1 — Ps)', the probability to produce a fast hadron 
from the background and not from hard process is Pb(1 — Ps)'i probability to produce a fast 
hadron in hard process and not in the background is ^^(l — p^), and PsPb is the probability 
to produce two fast hadrons - one in the background and one in the hard process. Since the 
last outcome contributes equally to Bpp and Spp we have 

Bpp = pb(1 - Ps/2), Spp = psil - Pb/2). (18) 

Since Spp is small, then to a very good approximation the solution of Eq. (fT8l) is pb = 
Bpp (1 + Spp/ {2 — B — S)) , ps = Spp (1 + Bpp/ {2 — B — S)). Hence pb is slightly larger than 
B, while for ps a relative correction is significantly larger. 

We can now calculate the probability that no hadrons will be produced in the inelastic 
collision of a nucleon with m nucleons of the nucleus: 



[l-B- S)^,ollisions = (1 - PbT{1 - Ps). (19) 

Using STAR data ioi S + B we find m = 2.8. It is easy to check that, due to ps <^ 1, 
this estimate of m is insensitive to the presence of two contributions to the multiplicity. 

The same picture allows one to estimate the value of S for dAu collisions. Qualitatively, 
we expect that S should drop as more hadrons are produced in soft collisions and the chance 
for the fastest hadron to be attributed to the recoiling jet becomes smaller. In the case of 
an inelastic collision of a nucleon with nucleons of the nucleus, the probability that in 
exactly m soft interactions a fast hadron would be produced, and that also a fast nucleon 
would be produced in a hard collisions is psC'^pBi)- — Pb)^~"^- For these events there is 
pa 1/(1 + m) chance that the fastest hadron would belong to the hard subprocess. Summing 
over m we obtain: 

m=N /^mf-\ „ \N—m^m 

c ST ^n{^-Pb) Pb fr,a\ 

ON collisions — Ps ■ 7 —r{ ■ l^UJ 

^0 + 1) 

Taking A^ ~ 3 we find S{dAu) ^0.1 which agrees well with the data. Thus we conclude 
that the increase of the associated soft multiplicity explains the reduction of S observed in 
the data without invoking any suppression of the recoil hadron production on the level of 
the hard subprocess. 

We have checked that accounting for the decrease of the soft multiplicity per Ncou leads 
to a small increase in our estimate (see also below). 

One can see from these equations that if the contribution of the central impact parameters 
{Ncoii ~ 13) were dominating in the vr*^ production like in the CGC models one would 
obtain {1 — B — S),S <^ 0.01 which is in a qualitative contradiction with the data. One 
would reach this conclusion even in the case of color transparency for the interaction of the 



"'^ Hereafter we are making an implicit assumption that one can neglect production of two hadrons from 
the soft or hard pp interactions within the experimental cuts. In the case of soft interactions this is justified 
both by small overall multiplicity and presence of short-range negative correlations in rapidities. In the case 
of hard process this is justified by a relatively small value of the pr of the trigger. Obviously one can improve 
this procedure by using information from the STAR experiment which is not available yet. 
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nucleoli involved in the hard interaction as the second nucleon would still experience ~ 6.5 
interactions. To compare predictions of the CGC models with data one should trigger for 
events at central impact parameters and look for suppression of recoil jets in such collisions 
using the method described in sectional 

Note also that a simple test of the relative importance of the central and peripheral 
mechanisms of the pion production is the ratio of the total hadron multiplicity in the events 
with the pion trigger and in the minimal bias events. We expect this ratio to be 

{-^trigger /-^min.bias) ~ (3/7.2) ~ 1/2. 

At the same time in the CGC model this ratio should be larger than one, since the relative 
contribution of the central impact parameters is enhanced as compared to the minimal bias 
sample by a factor A^^^ ~ 2.4 |13]. Unfortunately, information about this ratio was not 
released so far by the STAR collaboration. 

5 The distribution over the number of collisions 

In the previous section we calculated B and S for a fixed number of collisions. In a more 
realistic calculation we need to take into account distribution over the number of collisions. 
The important constraint here is that the suppression factor, RdAu, for inclusive 7r° produc- 
tion is RdAu ~ 0.3. This requires that at least nucleons with the impact parameter b > bmin 
satisfying condition 

J d%TA{b)e{b - bmin) = RdAu, (21) 

should contribute to the inclusive pion yield. Here T^(&) is the conventional optical density 
which is expressed through the nuclear matter density p^(r),/ d^rpA^r) = 1 as 

T4(6)= dzpAiVWT^). (22) 

J— oo 

Condition of Eq. (pTj) corresponds to 6 > 5fm for RdAu = 0.3. For collisions with 6 ~ 5fm 
the average number of collisions is already larger than 3. However, the presence of the more 
peripheral collisions still may lead to an average number of collisions close to 3. 

To simplify the discussion we will consider the case of pA scattering and later on correct 
for the presence of the second nucleon in the projectile. In the probabilistic/geometrical 
picture one can rewrite the inelastic cross section af^ as a sum of cross sections with exactly 
m inelastic interactions |25]: 

m=A ^1 

<^fn = E = • , , / d% (T^(6)ar )'"(1 - T^(&)^r )^-'"- (23) 

These partial cross sections satisfy the sum rule J2m=i ''^'^m = Aa^^. As a result, if the 
emission of the particles in each inelastic interaction is the same as in the NN collisions, all 
shadowing corrections are canceled reflecting the AGK cancelation [24]. 

To model distribution over the number of soft interactions, we need to introduce a sup- 
pression factor SF(b) which is a function of the nuclear density per unit area at a given 
optical density which is given by T(6). 
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Using this model we can calculate the A-dependence of the quantities measured by STAR 
taking into account the distribution over the number of the collisions: 



/-I R OX ^Z=l (1 - ^ - S) m collisions I 0'mib)SF{b)d b 

[l-B-b) = nnt^^o , (^4) 

^ in ^dAu 

where amib) is the integrand of Eq. fE^ and {1 — B — S)mcoiiisions is given by Eq. flTI?]) . 
Similarly, to calculate S{pAu) we can combine Eqs. [2Ul[2^ to find 

n 5^m,=l collisions 

Jam{b)SF{b)d^b 

^ in ^dAu 

For a numerical study we choose two models inspired by the energy loss estimate of 
section 2 for interaction near BDR and the regime of absorption deep in the BDR (cf. Eq. 

SF^'\b) = (1 + a,TAib))-\ SF^'\b) = (1 + a2TA(6))-^ (26) 



with the parameter ai = 2.5, a2 = 1.63 fixed by the condition j TA{b)SF''{b)d% = R^au = 
0.286 as measured by STAR for the higher pt correlation bin corresponding to averaging 
over 30 < i?^ < 55 GeV |7] . We found that the two models of suppression give very similar 
results for the observables measured experimentally with the second model for SF{b) giving 
slightly larger values of S and (1-B-S) since it yields a stronger suppression of scattering at 
the central impact parameters (the suppression factor is ~ 5.4 in the first model and ~ 20 in 
the second model). The numerical values are S=0.068 & .075; (1-B-S)=0.070 & 0.086. If we 
try to model the decrease of the multiplicity per wounded nucleon in line with the BRAHMS 
data we naturally find an increase of 5, (1 - 5 - S): S=0.085 & .090; (1-B-S)=0.11 & 0.12. 

However, in the actual experiment the dAu interaction was studied. In this case, the 
average number of collisions is about a factor 1.4-1.5 higher due to the interaction of the 
second nucleon. We can make a rough estimate of this effect by substituting collisions in 
pA scattering when the proton experiences N inelastic interactions by a superposition with 
equal probabilities of N and 2N inelastic collisions. Clearly, a more detailed modeling of 
dAu interactions is necessary - we will address it elsewhere. We find, when we account 
for the energy splitting S=0.067 & .072; (1-B-S)=0.066 & 0.079. These numbers should 
be compared with {1 — B — S) = 0.1, S = 0.093 ± 0.04 measured by STAR for the higher 
Pt bin which we analyze here. This suggests that, with inclusion of the second nucleon 
interaction in the Gribov-Glauber model, one gets a somewhat larger suppression of the jets 
than reported experimentally [f| and a smaller probability not to observe any hadrons than 
the one observed experimentally leaves room for effects due possible deviations from the 
geometrical picture. One effect of such kind which was suggested in [2S] is the presence of 
the color fluctuations the projectile nucleon. Selection of large xp in the nucleon may select 
fluctuations with smaller interaction cross section and lower the number of the interactions. 
However, the observables used in the analysis are not very sensitive to the distribution over 
the number of interactions as long as a peripheral model of tt^ production is used. 



^Use of two models allows us to test weak sensitivity of our conclusions to a choice of the specific model 
for dependence of suppression on the nuclear thickness. 

■^If some of the pions were due to a production mechanism without a recoil jet, S would decrease increasing 
discrepancy with the data. 
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To summarize, the results of our analysis of the data at higher px, the data are consistent 
with no suppression of the balancing hadron production for the trigger with (pt) ~ 1.3 GeV/ c 
0. Lack of the suppression of the pQCD mechanism for < xa >~ 0.01 which dominates in 
the correlation measurements of the STAR puts an upper limit on the x range where coherent 
effects may suppress the pQCD contribution. Since the analysis of pj[] find that the pQCD 
contribution is dominated by xa > 0.01, we can conclude that the main contribution both 
to inclusive and the correlated cross section originates from pQCD hard collisions at large 
impact parameters. 

The observation of the recoil jets in the pp case with a strength compatible with pQCD 
calculations suggests that the mechanism for pion production in the STAR kinematics is 
predominantly perturbative so that it is legitimate to discuss the propagation of a parton 
through the nucleus leading to pion production. To ensure a suppression of the pion yield 
at central impact parameters for the discussed kinematics one needs a mechanism which 
is related to the propagation of the projectile parton which is generating a pion in a hard 
interaction with the x ~ 0.01 parton. For example, the rate of suppression observed by 
BRAHMS would require fractional energy losses ~ 3% both in the initial and final state 
|llj . Similar losses would produce a suppression of the pion yield in STAR kinematics 
comparable with the inclusive data. Modeling performed above using Eg. (1261) indicates that 
for the central impact parameters the fractional energy losses should be at least a factor 
of 1.5 larger. Note here that such losses are sufficient only because the kinematics of the 
elementary process is close to the limit of the phase space. At the same time, this estimate 
assumes that fluctuations in the energy losses should not be large. For example, processes 
with energy losses comparable to the initial energy (like in the case of high energy electron 
propagation through the media) would not generate necessary suppression provided overall 
losses are of the order of few percent. Note also that the second jet in the STAR kinematics 
has much smaller longitudinal momentum and hence is far from the BDR. Therefore in the 
STAR kinematics one does not expect the suppression of the correlation with production of 
the second jet. However a strong suppression is expected for production of two balancing 
forward jets since both of them are interacting in the BDR. 

Hence the data are qualitatively consistent with the scenario described in the introduction 
that leading partons of the projectile (with x ~ 0.7) interact at central impact parameters 
with the small x nuclear gluon fields with the strength close to the BDR. 

6 Suggestions for future measurements to reveal onset 
of BDR 

We have seen that the quantities used in the STAR analysis involve an interplay of the hard 
and soft interactions. Here we want to suggest a few other observables which allow either to 
suppress this interplay or to optimize the sensitivity to the number of the collisions. 

First let us discuss another procedure for studies of the modification of the characteristics 
of the hard collisions which is significantly less sensitive to the properties of the soft inter- 

''in the case of the lower px trigger data set, our estimate for the dAu scattering gives {1 — B — S) = 
0.060 — 0.072 and S ~ .045 which is a bit above the reported value of = 0.020 ± .013. However application 
of hard scattering picture in any case rather problematic for pxitrigger) ^ 1 GeV/c. 
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actions. Let us consider the ratio of the double inclusive and single inclusive cross sections 
for production of a particle in forward and in central kinematics which are characterized by 
their rapidities and transverse momenta: 



where is the angle between— / and ptc- We can now introduce 

ARR{yf, \ptf\,yc, bid ,0) = RRivf, \ptf\,yc, IPtdA) - RRivj, \ptf\,yc, btd ,-0)- (28) 

Similar to the logic of the STAR analysis, we expect that only hard contributions to the 
central production depend on 0. Hence, in the case of inclusive quantities like ARR the 
soft interactions are cancelled, while (as we have seen above) this is not the case for the 
quantities considered in [7]. 

Consequently, the ratio of ARR'^^'^ and ARR^^ can be used to study how the pr, 0, t] 
dependences of the balancing jets depend on A (obviously one can consider the ratio of ARR 
integrated over all but one variable). The STAR analysis used a px > 0.5GeV/c cutoff to 
enhance the hard contribution. In our procedure one is likely to be able to use a smaller 
cutoff, or no cutoff at all. It appears that already current statistics of STAR would allow at 
least some of these measurements. Note here that the nuclear shadowing effects are more 
important for the positive rapidities of the recoil jets. Hence, a study of t] dependence in 
the kinematics studied by STAR could constrain the leading twist shadowing effects between 
0.005 < X < 0.02, albeit for rather large impact parameters where shadowing is smaller than 
on average. Note in passing that the current estimates of the suppression of the inclusive 
pion yield due to nuclear shadowing overestimate effect as they do not take into account 
that the process is dominated by the scattering at large impact parameters. 

It is worth emphasizing here, that for a large range of impact parameters, one is likely 
to be in the regime too close to BDR to apply the leading twist approximation for nuclear 
shadowing. At the same time the important feature of the leading twist nuclear shadowing is 
likely to hold, namely that Rg = Ga{x, Q'^)/AGn{x, Q"^) < 1 is achieved due to simultaneous 
interactions with 1/Rg nucleons, leading to an increase in hadron multiplicity at central 
rapidities and in the nuclear fragmentation region [27]. 

To study the dependence of pion production on the number of collisions, one needs 
to study the multiplicity distribution of the soft particle production at central rapidities. 
As a first step one would have to deconvolute the hard contribution which would be well 
constrained by a study of ARR (though this is actually a rather small correction, especially 
for large multiplicities). The tail of the distribution at large multiplicities will determine 
the relative contribution of the collisions with several nucleons - for the cuts of STAR the 
average multiplicity should grow ^ Ncou/'2,. 

This program would allow for a study of how the regime of large energy losses sets in as 
a function of the gluon / nucleon transverse density. Such a study would have important 
implications for LHC since in the large energy losses scenario, enhancement of the losses as 
compared to LT pQCD calculations is due to the proximity to the unitarity limit. Conse- 
quently, one would expect large energy losses for a much larger range of rapidities at LHC 
for the same parton virtualities. 




(27) 
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Two complementary methods to obtain information about the centrahty of dependence 
of the very forward pion production would be to use information from the zero degree 
calorimeters (ZDC) along the deuteron and gold beams. Measuring the number of neutron 
spectators produced in the fragmentation of the deuteron would be drastically different for 
the peripheral and central impact parameters scenarios - in the peripheral case one expects 
an increase of the spectator neutron multiplicity as compared to the minimal bias events 
as the neutron in peripheral interactions has a significant chance to survive (provided the 
pion was emitted in the interaction of the proton which occurs with 50% probability). We 
postpone a quantitative description of this significant effect for future publications. At 
the same time, a chance for a neutron to escape in the central collisions would be very 
small: exp(— cXm T{b ~ 0)) < 10-2 I. The ZDC measuring production of neutrons from 
the nucleus fragmentation would observe a smaller than minimal bias multiplicity for the 
peripheral scenario and a larger multiplicity in the central collisions scenario. An important 
advantage of these observables is that they are practically insensitive to the issue of the 
split of the energy between soft interactions. Hence, one can reduce uncertainties in the 
extraction procedure by employing information about the production of neutrons in generic 
dAu collisions, in particular in collisions with production of soft hadrons at central rapidities. 
It appears likely that such studies would substantially improve the determination of the T{b) 
dependence of the suppression factor. 

Future experiments at RHIC would allow one to separate large energy losses and leading 
twist nuclear shadowing. One would have to measure ARR as function of the rapidity. The 
shadowing effects would lead to drop of ARR'^^'^ (r]) / ARR^^ {i]) for forward rapidities where 
xa ~ 10^^ dominates. Also, if one would be able to decrease xn for fixed small xa, one 
would enhance the shadowing effects as compared to BDR effects. In this limit one would 
observe an increase in the associated multiplicity at the central rapidities since for the nuclear 
shadowing mechanism, central impact parameters give a large relative contribution. 

A color transparency effect would be manifested in a number of collisions with small 
Ncoii significantly larger than that given by the Glauber model. Obviously use of the proton 
beams would nicely complement studies with the deuteron beams as one would be able to 
compare triggers for centrality solely based on the interactions of one nucleon and on the 
interactions of two nucleons. 

After our first version of our analysis was completed PHENIX released the results of their 
analysis of the correlations [28] which are rather similar to the procedure we advocate. They 
study hadron correlations for smaller rapidities and higher pt than those studied by STAR 
which is far from the BDR and find no suppression of the correlations. 

7 Conclusions and open questions 

We have demonstrated that partons with large transverse momenta corresponding to rather 
large virtualities for which the BDR is reached should lose a substantial fraction of their 

^ Since the deuteron is weakly bound system, there is a significant tail in the wave function at distances > 
4fm making selection of pure central collisions sample very difficult. This problem would be greatly alleviated 
if one would study A-dependence of production of spectators with transverse momenta > 200 MeV/c which 
selects configurations with transverse separations < 2 fm. 
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energy. In the case of inclusive production of very forward pions this leads to the dominance 
of the scattering at peripheral impact parameters. For partons with xat > 0.5 the energy 
losses for central impact parameters should lead to suppression of the inclusive yield at least 
by a factor of five which corresponds to energy losses > 10%. As a result BDR leads to an 
extension to a wider pt range of the pattern of the strong suppression of the leading hadron 
production at small pt observed in the central pA collisions. 

With increase of energy from RHIC to LHC energy losses at large should strongly 
increase, while substantial losses > 10% should persist for rapidities \y\ > 2. It appears that 
this should lead to increase of the densities in the central collisions as compared to the current 
estimates. It will also lead to suppression of the production of the recoil jets at the rapidity 
intervals where no suppression is present at RHIC. In the forward direction we expect a 
significantly larger suppression than already large suppression found in ^20j where fractional 
energy losses were neglected. Fractional energy losses result in modification of the form of the 
QCD factorization theorem at LHC energies. In particular they lead to suppression of Higgs 
/ SUSY particle production in pp scattering by at least [(1 — x/(l + e))/(l — x)]^^ . Here 
X > vtijil^fs are the light-cone fractions carried by initial gluons which initiate production 
of the Higgs particle with accompanying bremstrahlung, and e > .1 is the fractional energy 
loss. This corresponds to a suppression > 10% for mu = 140 GeV. 

Further studies of the proton/deuteron - nucleus interactions at the central impact param- 
eters at RHIC and future experiments at LHC would provide important constrains on this 
important ingredient of high energy dynamics. Similar effects will be present in the central 
pp collisions at LHC. They would amplify the correlations between the hadron production 
in the fragmentation and central regions discussed in Ref. [8]. 

We would like to thank Yu.Dokshitzer, S. Heppelmann, D. Kharzeev, L. McLerran, 
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